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The polarized Gowdy T 3 vacuum spacetimes are characterized, modulo gauge, 
by a "point particle" degree of freedom and a function ip that satisfies a linear 
field equation and a non-linear constraint. The quantum Gowdy model has been 
defined by using a representation for tp on a Fock space T . Using this quantum 
model, it has recently been shown that the dynamical evolution determined by 
the linear field equation for tp is not unitarily implemented on T . In this paper: 
(1) We derive the classical and quantum model using the "covariant phase space" 
formalism. (2) We show that time evolution is not unitarily implemented even on 
the physical Hilbert space of states TL C T defined by the quantum constraint. (3) 
We show that the spatially smeared canonical coordinates and momenta as well as 
the time-dependent Hamiltonian for tp are well-defined, self-adjoint operators for 
all time, admitting the usual probability interpretation despite the lack of unitary 
dynamics. 



1. Introduction 



Over the past 30 years, spacetimes admitting two commuting Killing vector fields 
have been studied repeatedly as "midi-super space" models for canonical quantum gravity; 
see, for example, references [1-11]. These models admit an infinite number of degrees of 
freedom — they are field theories — and as such they are more sophisticated than the 
"mini-superspace" models, which are mechanical models with a finite number of degrees 
of freedom. In particular, the quantized midi-superspace models can bring into play in- 
trinsically quantum field theoretic features which have no analogs in quantum mechanics. 
One of these features - the failure of time evolution to be unitarily implemented - is the 
impetus for the present paper. 

Of the midi-superspace models the Gowdy class [12] is interesting since it defines 
an inhomogeneous cosmology including a big bang or a big crunch. It was first studied 
as a model of quantum gravity by Misner [2] and Berger [3], who explored a variety of 
approaches to defining the quantum theory and extracting physics from the polarized 
Gowdy T 3 model. This model arises by assuming spacetime is not flat, that it has the 
topology R x T 3 , and that there is an Abelian 2 parameter isometry group with spacelike 
orbits T 2 C T 3 generated by a pair of commuting, hypersurface-orthogonal Killing vector 
fields. In this setting, the vacuum Einstein equations imply that, modulo gauge and a 
"point particle" degree of freedom, the classical dynamics of this model is governed by 
a function ip that satisfies a linear field equation and a non-linear constraint. The field 
equation is equivalent to that of a one-dimensional symmetry reduction of a massless 
free scalar field propagating on a flat spacetime (M,g). The spacetime (M,g) can be 
viewed as the causal region of (a compactification of) a three-dimensional version of Misner 
spacetime.* The Gowdy time foliation equips the (M,g) with a foliation by expanding (or 
contracting) spatial sections, along which the time evolution of <p is given by the linear field 
equation. The constraint requires the total momentum of the scalar field <p to vanish. One 
of the quantizations of <p studied by Berger, and subsequently studied by Husain [4], Pierri 
[10], and Corichi et al [11], is equivalent to the restriction of the standard quantum theory 
of a massless free scalar field on three-dimensional Minkowski space to the compactified 
Misner spacetime (M,g). This defines a Fock space representation for ip. The total field 
momentum can be defined in this representation and physical states are eigenstates of it 
with zero eigenvalue. Thus the polarized Gowdy T model is defined as a constrained 
quantum field theory. 

In reference [11], it was observed that the time evolution of <p, as defined by its linear 
field equation, cannot be implemented as a unitary transformation on the Fock space 
described above. This sort of phenomenon, which is not unexpected in quantum field 



* For a discussion of Misner spacetime, see [13]. 
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theory [14,15], has been seen in other, related settings [16,17]. The lack of unitary dynamics 
leads the authors of [11] to conclude that the quantized model under consideration is not 
physically viable. In this paper we will extend the results of [11] in 3 significant ways, 
which will be described in the following paragraphs. 

First, the classical analysis we develop in §2 as the underpinning for the quantum 
theory (§3) is quite different from the standard Dirac techniques utilized in [10,11] since 
we formulate the model using the "covariant phase space" formalism [18]. This feature of 
our work is perhaps of some intrinsic interest since it represents a non-trivial application 
of that formalism. But the principal utility of the covariant phase space approach to the 
Gowdy model is that it allows for an independent, relatively simple construction of the 
"deparametrized" dynamical system, which is not entirely straightforward in the Dirac- 
type of approach featuring in [10,11] owing to the presence of the point particle degrees of 
freedom which are subsequently mixed in with the time variable during the deparametriza- 
tion process. The covariant phase space formalism allows one to work directly with the 
classical spacetime and this makes it very easy to keep track of the field degrees of freedom, 
the point particle degrees of freedom and the choice of time that is used in the model. 

Second, we show in §4 that the proof of the failure of time evolution to be unitarily 
implemented on the initial, "auxiliary" Fock space can be extended to apply to the physical 
Hilbert space of states defined by the momentum constraint. A priori, it is possible that 
the dynamical evolution is unitarily implemented on the physical Hilbert space but not on 
the auxiliary Hilbert space. In this scenario the polarized Gowdy T 3 model would have 
unitary dynamics, physically speaking, with the non-unitarity found in [11] just being a 
technical complication. We show, however, that this is not the case. 

Third, we show in §5-6 that, despite the failure of time evolution to be unitarily im- 
plemented, a number of basic operators have entirely satisfactory behavior. In particular, 
the canonical coordinates and momentum of the quantum field <yC, and all their derivatives 
(when spatially smeared with smooth functions) are at each time well-defined self-adjoint 
operators with continuous spectrum on the real line. Therefore, if not for the existence of 
the momentum constraint, which restricts the class of physical operators, the canonical field 
operators would represent observables with a perfectly acceptable physical interpretation 
in the Heisenberg picture. Operators representing physical observables can be obtained 
from the canonical field operators by projection into the physical Hilbert space [10], and 
the self-adjointness of the field operators implies that these operators are again physically 
acceptable, despite the lack of unitary dynamics on the physical Hilbert space. Moreover, 
we show that the one-parameter family of Hamiltonians for the deparametrized Gowdy 
model can be defined as self-adjoint operators both on the auxiliary Fock space and on the 
physical Hilbert space, despite the fact that the dynamical evolution they generate is not 
unitarily implementable. 
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The results summarized in the previous paragraph indicate that it may be possible to 
considerably soften the conclusions reached in [11] concerning the physical viability of the 
quantum Gowdy model. We discuss the situation in §7. 



2. Polarized Gowdy model: classical theory 

Fix coordinates (t, x, y, z) on M = R+ x T 3 , where t > and < x, y, z < 2tt. The 
polarized Gowdy T 3 metrics are defined by [12,19] 

g = I j e (7-¥>)(_<ft ^ dt + dx <g> dx) + t^e'^dy <g> dy + e^dz <g> dz} , (2.1) 

where / > is a constant and both 7 = ^(t, x) and ip = (p(t, x) are periodic functions of x 
with period 27r. The space of Gowdy metrics is thus parametrized by (l,^,ip). 

The vacuum Einstein equations, when restricted to (2.1), are equivalent to 

-<P,tt ~ j<P,t + ¥,xx = 0, (2.2) 



7,t = g + V 2 x) = l,x = t(p,t<P,x = n. (2.3) 
All smooth solutions to (2.2) are of the form 



1 1 00 

= + pint) + -= E («niyo(|n|*)e^ + a*^ *(|n|t)e-^), (2.4) 

V27il 2V21 v / 



n— — 00 



where ifg i s a Hankel function of the second kind* and the sequence {a n } is rapidly 
decreasing, i.e., its elements approach zero faster than the reciprocal of any polynomial in 
nasn-> ±00. Given (2.4), the solutions to (2.3) can be expressed in the form 

7(i,x)=70+ / dt'£(t',x)+ / dx'U(t ,x'), (2.5) 
J to Jxq 

where 70 is a constant. The metric function 7 is thus parametrized by 70 and cp. (The 
quantities (to, %o) can be fixed arbitrarily; different choices of them merely redefine 79.) 7 
is periodic in x if and only if 

C = l dxU(x)= ri\a n \ 2 = 0, (2.6) 

^° n=-oo 

* We could equally well use Hankel functions of the first kind, which is the choice made in 
[10,11]. The choice used here has the feature that its notion of positive frequency agrees 
with that defined by the timelike Killing vector field of Minkowski space when the field 
ip is interpreted as a function on the interior of the future light cone. This is the same 
convention used in the "Schmidt model" [20] . No results in this paper depend upon which 
type of Hankel function is used. 



which can be viewed as the sole remnant of the Hamiltonian and momentum constraints. 
Note that C is independent of t, so that the constraint (2.6) need only be imposed at a 
single value of t. 

According to (2.4)-(2.6), the set of smooth solutions to the Einstein equations for the 
polarized Gowdy metrics can be identified with the set 

r = {(l, 10 ,q,p,a n ,a* n ), n = ±1, ±2, . . . \C = 0}. (2.7) 

The set V has a pre-symplectic structure naturally defined by the Einstein-Hilbert action as 
follows.* With a convenient normalization, the Einstein-Hilbert action takes the following, 
remarkably simple form when restricted to the Gowdy metrics (2.1): 




Of course, one cannot obtain all the field equations (2.2)-(2.3) as critical points of (2.8) 
since gauge fixing conditions have been used to define (2.1) .f Nevertheless, as we shall see, 
this action is suitable for defining a pre-symplectic structure on r.f Following the general 
prescription of [18], vary the action (2.8) and extract the resulting boundary term, which 
defines a (t-dependent) 1-form on the space of dynamical variables (Z, 7, ip). Pull back this 
1-form to T to get the pre-symplectic potential. At a point (Z,7o, <p) G T, the value of this 



* Using the results of [21] it can be shown that the polarized Gowdy isometry group sat- 
isfies the "Principle of Symmetric Criticality" . This implies that any generally covariant 
local variational principle for a spacetime metric, when restricted to the set of all met- 
rics admitting the polarized Gowdy isometry group, yields the correct set of reduced field 
equations. 

f One can get only the p field equations from varying p in (2.8). The equations obtained 
from varying 7 are trivial since 7 only appears as a divergence. The equation coming from 
varying Z is non-trivial, but is automatically satisfied when p and 7 satisfy the equations 
(2.2)-(2.3). 

f From the point of view of a traditional Dirac type of Hamiltonian analysis, one can interpret 
this action and its pre-symplectic structure as corresponding to the result of a partial 
gauge fixing and deparametrization, with all of the first class constraints except (2.6) 
being rendered second class, followed by the use of Dirac brackets. 
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1-form on a tangent vector (51, Sjq, Sip) is given by 

"27T 



8(61,670, M = / dxl (tS^ t -t6(p t - 6~/ + t(p t S<p) 
Jo 

r2ir f rX rt 

= dxll[-6'yo + t6£ -tS(p tt - dx' 6U(t , x) - j dt'6£(t' ,x) + t<p >t 6(p 

JO I Jx Jto 

= I ^ - 2tt57o + t6E(t) - V2^6 (j^J - J dx J dx'5U(t , x) - dt'6E(t' 



where 



r2n 

+ J dxt<pj6<p 

= 16 ^-2tt7 + tE(t) - )JyP ~ J dx j dx'U(t ,x') - dt'E(t') 
p2tt 

+ dxl tip t5<p 
Jo 

(2.9) 

p2tt 

E(t) = / dx£(t,x). (2.10) 
Jo 

One point to keep in mind here is that the fields <p and their variations 6<p are subject 
to the field equations and constraints. The pre-symplectic two-form is obtained from the 
symplectic potential by exterior differentiation. At the point (l,jQ,ip) e T, the two-form 
evaluated on a pair of vectors, (SI, Sjq, Sip) and (61, 5^q, Sep), is given by 

u(6l, S-fo, Sip; 61, 6-fo, Sep) = SIS£ - SIS£ + J dx (sP^ty - SP^S^ , (2.11) 

where 

:= Vlp, P^ := tVlpj, (2.12) 

and 



/2-Ti- /• 27r r x ft 
\ —p- / dx / dx'U(t ,x')- / dt'E(t') 
V I Jo Jx a Jt 

1 1 f Z7T 
+ 2lJ ^^(^W'^O- 



£ : = -2n 10 + tE(t)- X l—p- 



(2.13) 



From the field equations it follows that £ (like I and 70) is a constant of motion: 

ft = °- < 2 - 14 > 

From the field equations and their linearization it follows that 

^u(6l, Sjq, Sip; 51, SjO, Sip) = 0, (2.15) 



so the pre-symplectic structure is defined independently of the value of t. This is guaranteed 
by the nature of its construction [18].* 

The form to is degenerate precisely because ip is subject to the constraint (2.6). The 
degeneracy directions for u> are given by 

61 = = 5£, 5<p = <p, x , (2.16) 

which implies, in particular, 

5 10 = n x (to,x ), (2.17) 

so that, 

uj(81, ^70, <fy; 0, n,a:(*o> xo), (f tX ) = 0. (2.18) 

The vector field on V defined by (2.16) generates a 1-parameter pre-symplectic group of 
transformations on T: 

<p(t t x)^<p(t,x + \), l^l, (2-19) 

We can interpret the pre-symplectic structure on V as follows. Consider the uncon- 
strained space r = {(/, 70, q,p, a n , a^), n = ±1, ±2, . . .}. Define a non-degenerate 2-form 
Co on this space using (2.11), but without imposing C = = SC. It is straightforward to 
verify that Co is conserved by virtue of the field equations and their linearization, so that 
this symplectic form is t-independent. From Co we can read off a canonical chart for T; the 
canonical pairs are ((f), P^) and (£, /). Since / > it is convenient to define new canonical 
variables (Q, P) G R 2 by 

Q = ]nl, P=-l£. (2.20) 

We extend the one parameter group (2.19) and its infinitesimal generator (2.16) to V in 
the obvious way. The pre-symplectic space of solutions (I\a;) can then be re-constructed 
from (r,a>) by (i) imposing the constraint, 

p2tt 

C = ^ P^, x = 0, (2.21) 

thus defining V C T, and (ii) pulling back the symplectic 2-form Co to T, thus defining uj. 

To summarize, the space of solutions to the vacuum Einstein equations for the Gowdy 
metrics (2.1) consists of a "point particle" degree of freedom with canonically conjugate 
variables (Q, P), and a field degree of freedom described by (p, or equivalently the canonical 
variables (0, P^), or equivalently the variables (q,p, a n , a^) n = ±1, ±2, . . ., subject to the 

* In fact, a similar computation shows that -^9(51, Sjq, dip) = 0. 9 is conserved because of 
a special feature of the Einstein-Hilbert action: it vanishes when the field equations are 
satisfied. 
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constraint (2.6), or (2.21). This characterization of the vacuum Gowdy spacetimes is 
essentially the same as obtained from the Hamiltonian methods of [7,10,11]. The only 
significant differences are as follows. First, by working directly with the space of solutions 
to the field equations, the metric variable 7 can be defined in terms of tp using an integral 
over time (and space) rather than the purely spatial integral that is used in the Hamiltonian 
formalism. Second, the "point particle" degrees of freedom (Q, P) have been defined as 
constants of the motion and have been clearly disentangled from the time variable. There 
is no "q-number" aspect to the time, such as arises in [10,11].* 

Dynamical evolution from t = t\ to t = t 2 can be viewed as a pre-symplectic map 
'It 1 t 2 :T — > r (see, e.g., [17]) . The map takes a given solution (l,"fQ,(p) of the Einstein 
equations to a solution (£,70, <£) whose Cauchy data at t = t\ are the Cauchy data for 
(I, 70, if) at t = t 2 - In particular: 

m) = (to) 

h<fi,t(tl) = *2<P,t(*2)- 

Since the general solution to the field equations is known explicitly, it is straightforward 
to construct (l,%,0) = T^j^l, 70, <p). We have 

70 = 70, 1 = 1 ^ Q = Q, P = P, (2-23) 
0(t,x) = -=(q + plnt) + —= (a n H (\n\t)e inx + a* n H* (\n\t) e - inx ), (2.24) 



2V21 n ±_ x 

n#0 



where 



and 



with 



and 



t7T . . 

a n = -j\n\ 



t 2 

q = q + pln(— ), p = p, (2.25) 
h 



a n = OL n a n + f3 n at n , (2.26) 
hHKlnltjHoQnfo) - t 2 H^\n\t 1 )H 1 (\n\t 2 )] (2.27) 



a m \ 1 
Pn = -£\n\ 



hHtdn^HZdnlti) - ^H^n^Ht^) \ . (2.28) 
The same equations, (2.23)-(2.28) define a symplectic map T: V — > V. 



* These authors rescale the time variable with a dynamical variable constructed from the 
point particle degree of freedom so as to simplify the form of the field equations for the 
field (p. However, this complicates the dynamical behavior of the point particle degrees of 
freedom. By working with constants of motion, these complications are avoided. 
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3. Polarized Gowdy model: quantum theory 



The quantization of the polarized Gowdy model used in [3,4,10,11] can be understood 
in the present formulation of the model as follows. Define a Hilbert space H by 

H = L 2 (R 2 )®F, (3.1) 

where JF is the symmetric Fock space built from the Hilbert space of square-summable 
complex sequences, ( n , n = ±1, ±2, ... , 

J2\(n\ 2 <oc. (3.2) 

Any $ G J can be represented as an infinite sequence of complex sequences* 

* = (V'0»V'mi»V'mim2>---»V'mi-mfc>---)> ( 3 - 3 ) 

where ipQ G C, 
and 

oo 

l^0| 2 + $^ ^ l^mr-mj 2 < OO. (3.5) 
fc=l mi-rajj^O 

The canonical pairs (Q,P) and (g,p) are represented as identity operators on JF and are 
represented on (a dense domain in) L 2 (R 2 ) exactly as one would canonical coordinates 
and momenta for a particle moving in two dimensions, e.g., 

2 2 1 1 

ip = ip(x, y) e L (R ), qip = xip, pip = -d x ip, Qip = yip, Pip = -d y ip. (3.6) 

% % 

(Other equivalent representations are, of course, possible.) The remaining degrees of free- 
dom in the field if are represented as identity operators on L 2 (R 2 ) and are represented 
on the Fock space T as follows. Using the representation (3.3) for $ 6 f, we define 
annihilation and creation operators for each / ^ by 

ai^ = (iphV2iP lmi ,V3ip lmim2 ,...) (3.7) 
tfV = (0, ipoS mih V28^ mi ip m2 y V35^ mi ip m2m3 y ...). (3.8) 

These operators (on their common domain) satisfy 

a n = ( a nV, [«n,4] = &nml- (3.9) 

* Here ipQ is the "vacuum amplitude" , ip mi is the amplitude for "1-particle with momentum 
mi", etc. For convenience, we represent the entire sequence {ipk, k = ±1, ±2, . . .} simply 
by the symbol ip^. 
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The quantum field ip is denned as an operator- valued distribution on R+ x S 1 using the 
operator representation of (q,p, a n , a*J in the expansion (2.4). This quantization just 
described satisfies the prescription 

{Poisson Bracket} <-> - {commutator} 

% 

for the canonical coordinates on f, where the Poisson algebra is defined by the symplectic 
form uj. 

The representation of (a n , a*J on T just described can be viewed as the outcome of a 
general procedure, in which the Fock space representation is defined once the appropriate 
"one particle Hilbert space" is extracted from the space of solutions of the equation (2.2) 
(see, e.g., [22]). Let us sketch this construction since some of its ingredients will be useful 
in what follows. Modulo the zero frequency modes (q,p), the space of solutions to (2.2) 
can be identified with the space S of rapidly decreasing sequences of complex numbers 
p = {p n }, n = ±1, ±2, . . .. The symplectic form u can be pulled back to give a symplectic 
form O on S. Introduce a scalar product p: S x S — > R by 

Mp> a ) = \ ^2(Pn a n + (TnPn)- (3-10) 

Following the prescription found in reference [22] , this scalar product can be used to define 
a Hilbert space H of square summable sequences of complex numbers (denoted as before) 
such that the Hilbert space scalar product is given by* 

i 

(p,a) = p,(p,a) - -Q(p,a). (3.11) 

H is the one particle Hilbert space out of which the symmetric Fock space JF is constructed 
via tensor products and direct sums, as usual. 

The relation between this way of describing the Fock representation and that of refer- 
ences [10,11] is as follows. In references [10,11] the one particle Hilbert space is extracted 
from the space of solutions S of (2.2) by defining a suitable complex structure J:S — > S. 
This complex structure can be obtained from the scalar product p and symplectic structure 
O by "raising an index" on Q using p so that 

n(p,a) = 2p(p,Ja). (3.12) 

As noted above, this quantization is based upon the canonical commutation relations 
associated to the symplectic space (r,<2>), which is the space of solutions to all field equa- 
tions except for the constraint (2.6). The usual strategy for imposing the constraint is to 

* The bilinear forms p and O on S are extended to H by complex linearity and continuity. 
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represent C as a suitable operator on Tt and then to define the Hilbert space of physical 
states, T~t p hys, as the kernel of C. To do this we proceed in a slightly roundabout fashion 
which is technically more convenient and which provides a simple example of a general 
strategy we shall use again when discussing the Hamiltonian. Essentially, we shall define 
C as the generator of the unitary transformation that implements the classical transfor- 
mation (2.19). The physical Hilbert space is then defined as the set of vectors which are 
invariant under the unitary group. 

The transformation (2.19), extended in the obvious way from V to T, corresponds to 
the symplectic transformation 

(q,p)^(q,p), (Q,P)^(Q,P), (a n ,a* n )^(e mX a n , e- mX a* n ). (3.13) 

It is easy to check that this 1-parameter group is strongly continuous* in A relative to the 
norm 

\\p\\ 2 = p(p, p) (3-14) 

defined on H. It then follows from the results of [15] that the symplectic group (3.13) 
is implemented in the sense that there exists a (strongly) continuous group of unitary 
transformations U(\): T — > T - uniquely determined up to a phase for each A - such that 

U\X)a n U(X) = e mX a n . (3.15) 

In fact, since the Fock representation is the GNS representation of a state canonically 
built from \x [22], unitary implementability follows directly from the observation that the 
symplectic transformation preserves the scalar product p [23]. Choosing the phases so 
that the Fock vacuum state is invariant under the unitary group, we have, using the 
representation (3.3), 

[/(A)* = (^o, e- 8fflI V mi ,e-' (mi+m2) V mira2 , • • • , e-^ +m * + --- +m ^ mim2 ... mk , ...). 

(3.16) 

The quantized constraint (again denoted C) is defined as 



.dU{\) 
L = i 



(3.17) 

A=0 



dX 

Using the representation (3.3), C is given (on an appropriate dense domain) by 

C^> = (0,mi^ TOl , (mi + m 2 )'i/ , m 1 m 2 , • • • , {mi + m 2 H h m k )ip mi m 2 ...m k , . . .). (3.18) 



* The transformation group p — > p(X) is strongly continuous if and only if lim^^^ o ||p(A) — 



p(Ao)|r = 0, VA GR. 
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We extend the definition of U(X) and C to H by defining these operators to be the identity 
on the L 2 (R 2 ) factor. We define the physical space 'Hphys C H as the set of vectors 
invariant under the unitary group U(X). 

Kphys = G H\U{\)V = VA G R}. (3.19) 

It follows that \1/ G Hphyg if and only if its Fock components (3.3) satisfy 

ip mi = 



^mim 2 = <J(mi + m2)xmi, 
^ mi m 2 ---m k = S(mi + ra 2 H h m^Xmimj-mn-n 



(3.20) 



where the sequences {V'mi > ^ ) m 1 ra 2 % rt Pm 1 m 2 ---m k } are each square-summable. This is formally 
equivalent to the usual definition of 'Hphys m which C is defined by normal-ordering and 
the physical states are annihilated by C. 

4. Non-unitary dynamics on Ti-phys 

The principal finding of [11] is that the Gowdy model time evolution is not unitarily 
implemented on the auxiliary Hilbert space 7i. In our formulation of the model this result 
arises as follows. 

Time evolution is defined in the Heisenberg picture by the (pre-)symplectic transfor- 
mation (2.23)-(2.28). Time evolution is unitarily implementable on 7i if and only if there 
exists a unitary transformation U = U{t\, ty): H — > H such that 

U^QU = Q, U^PU = P, (4.1) 

tfqU = q + pln(—), tfpU = p, (4.2) 
h 

U^a n U = a n a n + f3 n aL n , U^a^U = a* n ah + Pn a -n- (4.3) 

There is no obstacle to satisfying (4.1) and (4.2), but (4.3) is possible if and only if the 
sequence {/3 n } is square-summable [24] (also see [22] and references therein). From the 
large-argument asymptotic expansions of the Hankel functions appearing in [3 n it follows 
that f3 n is not square-summable, 

\Pn\ 2 = + O(^) =► £ l^l 2 - oo, (4.4) 
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so U, as defined above, cannot exist. 

Strictly speaking, this does not establish that the time evolution fails to be unitarily 
implemented in the Gowdy model because the Hilbert space Ti is merely an auxiliary device 
used to construct the physical Hilbert space 7~t p hys- Time evolution need only be unitarily 
implemented on Hp^ys'i a priori it is possible that a unitary U : Hphys ~^ ^phys exists, but 
has no appropriate extension to all of H. As it turns out, this scenario does not occur and 
time evolution cannot be unitarily implemented on Hphyg either. So this loophole in the 
argument for non-unitarity given in [11] can be closed. We demonstrate this as follows. 

We first note that the product of any vector x £ L 2 (R?) and the "vacuum state" in 

* = X® (1,0,0,0,...), (4.5) 

is a "physical state" , that is, 

U(\)* = *0, (4-6) 

so that £ 1~Lphys- Next we note that a_ k a k , k ^ 0, defines a linear operator on (a dense 
domain in) 'Hphys- We define the time evolution of this operator to be that induced by the 
Bogoliubov transformation (2.26)-(2.28): 

«-fc«fc -»• (&_ k a_ k + 0_ k a\){a k a k + Pk a _ k )- ( 4 -7) 

It is easy to check that the right hand side of (4.7) defines a linear operator on (a dense 
domain in) T~tphy S - We now show that the square-summability condition on {/3 n } is again 
necessary for the transformation (4.7) to be unitarily implementable on 'Hp\ l y S - To this 
end, we suppose that there is a family of unitary transformations 

U(ti, t 2 ) : 'Hphys — > Hphy S (4-8) 

that is continuous for all t\ , t 2 > and satisfies 

U(t, t) = I, U(ti,t 2 )U(t 2 , h) = U(h, t 3 ), (4.9) 

U Jf (ti,t 2 )a_ k a k U(t 1 ,t 2 ) = (a_ k a_ k + P- k a\){a k a k + /3 k a}_ k ). (4.10) 
The condition (4.10) implies that 

(a_ k a_ k + (3_ k a\){a k a k + P k al k )u\t 1 ,t 2 )^o = 0. (4.11) 

From (4.6) and (3.20) we have that U^q takes the form 

W(ti,t 2 )^o = x® (40,<5(mi +m 2 )ip mi ,S(mi + m 2 + ms)ip mi m 2 , . . .), (4.12) 

for square-summable ip mi , / 4>m 1 m 2 i ■ ■ ■■ The vacuum component of (4.11) implies 

Vfc = -^Vo- (4-13) 
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From (4.13) and (3.5), a necessary condition for U^{t\, £2)^0 ~~ an d hence U(t\,t2) - to be 
denned is square summability of {V'fc}- Because the sequence is bounded from below 
away from zero and also bounded from above, square-summability of {ipk} ls equivalent 
to square summability of assuming that ipQ 7^ 0. Because U(ti,t2) is continuous in 

each argument, it follows from (4.9) that ipQ 7^ 0. Therefore, unitary implementability on 
'Hphys requires square summability of which does not hold, (4.4). 



5. Self-adjointness of canonical field operators 

Here is is shown that, despite the failure of time evolution to be unitarily implemented 
on 7i or Hphy S , the canonical coordinates and momenta for the field tp are well-defined, 
self-adjoint operators for alH > and admit the usual probability interpretation. 

The canonical field operators (</>, Pa) associated to a time t are formally defined as 
distributions on S 1 via (cf. (2.4) and (2.12)) 

1 1 00 

(j) ( x )= ( q + plnt) + - 7 =Yl (mn\t)e inx a n + HZ(\n\t)e- inx aV), (5.1) 

* * n— — 00 

1 ,00 

W = V^ P ~ 271 ^ H (^(N*) 6 "^ + Hf(\n\t)e- inx aV) • (5.2) 



n— — 00 
n#0 



These distributions formally define canonical field operators (</>(/), P^ig)) associated to 
any smooth real- valued functions / and 5 on S 1 : 



<!>(/) = h(q + pint) + Y- (H (\n\t)f- n a n + H$(\n\t)f n aV), (5.3) 



n— — 00 
n#0 



P (j>{g)=gOP l n ' (■ ff l(l n l*)fl , -nan + -ffl(l«l*)fl'nOn) , ( 5 - 4 ) 

where 



n— — 00 
n#0 



OO OO 



-I — 1 ^ 



27T ^— ' V27T 

Ti= — oo v n= — 00 

Note that the sequences {f n } and {g n } are rapidly decreasing. 

To make these formal definitions precise we consider the dense subspace 

H = Dxf fl cK, 

where V C L 2 (R 2 ) is a common, invariant dense domain for the operators (q,p,Q,P) 
(e.g., in the representation (3.6) V could be chosen to be smooth functions with compact 
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support) and J-'q C T is the dense set consisting of vectors with a finite number of non-zero 
components when represented according to (3.3) (i.e., Fock states with "a finite number of 
particles"). It follows that Ho is a dense, invariant domain for the canonical field operators. 
Moreover, it is easily checked that these operators are Hermitian (i.e., symmetric) on Hq. 

Because q and p are already defined as self-adjoint operators which commute with a n 
and an, in order to show that </>(/) and P^g) are self-adjoint it is sufficient to restrict 
attention to the field operators modulo the constant modes, which amounts to using test 
functions / and g in (5.3) and (5.4) with /o = = 90- We then proceed by using the 
approach described in §X.7 of [25]: prove that J-'q is an analytic domain for <p(f) and P(h(g) 
(without the constant modes). Nelson's analytic vector theorem then implies essential self- 
adjointness of the canonical field operators and hence a unique self-adjoint extension. The 
details follow. 

Begin with a Fock state with exactly N "particles" : 

^ = (0,0,...,^ mi ... mw ,0,0,...) GJ^o- 
With t fixed but arbitrary, define 

a(fH ) = f-nH (\n\t)a n , a\fH Q ) = £ f n H* k (\n\t)a{ 

It is straightforward to verify the inequalities 

\HfH fn < ^^±^.\\fHo\m\\, \\a\fH y*\\ < ^±^\\fHo\m\\, 

(5.6) 

where 

||/iy || 2 = ^|/n^0(Ht)| 2 , |H| 2 = l^mr-mf. (5.7) 

n^O mi--m„^0 

The first sum in (5.7) converges provided {f n } is square-summable, which it is since /: S 1 — > 
R is smooth. The second sum in (5.7) converges since {i)m 1 --m n } must be square-summable 
if \1> is to be in H. The inequalities (5.6) guarantee that 

\\Hfrn < ^^^\\f H o\\ ihi, (5.8) 

By definition, \1/ is an analytic vector for <p(f) when 

oo i 

;=o L 

From (5.8) it follows that 

ll^*H 7[ ^ ^^npWoll IMI, (5.io) 
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from which it follows that (5.9) is satisfied. This result easily generalizes to superpositions 
of vectors with different (but finite) numbers of "particles" and thence to all of JF . 

The preceding paragraph shows that, with respect to </>(/), Hq is a dense, invariant, 
analytic subspace of H. Therefore </>(/) is essentially self-adjoint on this domain and has 
a unique self-adjoint extension [25]. In a similar fashion one can show that P^(g) can be 
defined, for all time t > 0, as a self-adjoint operator for any smooth function g. The only 
change needed in the argument given above is the replacement 

#0(N*)/n^*M|#l(Mt)2n. (5.11) 

Indeed, as long as the spatial smearing functions are smooth, all spacetime derivatives of 
ip can, for each t > 0, be defined as self-adjoint operator-valued distributions on S . 

The spectrum of each of the canonical field operators and their derivatives is the whole 
real line, for all t > 0. This can be seen, for example, by comparing expressions such as (5.3) 
and (5.4) with their counterparts for a free field on a flat spacetime in inertial coordinates, 
which have continuous spectra for all choices of the test functions. The Hankel functions 
appearing in (5.3) and (5.4) can then be viewed as simply redefining the test functions. 

The spectral theorem [26] guarantees that for each self-adjoint operator A on a Hilbert 
space TC there is a unique-projection- valued measure cr(fi) associated to any measurable 
set O C R such that 

A= I \da{\), (5.12) 

Jr 

Given a state represented by the unit vector \I> e H, the probability Va{Q) that the 
observable represented by A is found to take the value in the set O C R is given by the 
expectation value 

V A (n) = (y,<r(m)- (5-13) 

Temporarily ignoring the constraint and working on H, this result can be applied to the 
(spatially smeared) canonical field operators at each time t > 0, whence the usual prob- 
ability interpretation can be implemented. Note, in particular, that probabilities always 
add up to unity because of the general property 

<r(R) = /, (5.14) 

where / is the identity on 7i. 

Of course, H is not the physical Hilbert space Hphyg, nor do the canonical field op- 
erators represent observables since they are not linear operators on V, v hy S . The point is, 
however, that the failure of unitary implementability of dynamics on Tt does not destroy 
the physical viability of the field operators (ignoring the constraint) . Taking account of the 
constraint is technically more complicated, but does not alter this conclusion. Self-adjoint 
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operators on H-phys can be defined by composing polynomials in the (spatially smeared) 
canonical field operators at each time t > with the projection operator into Hphys [10]- 
From the spectral theorem, these operators represent observables with the usual probabil- 
ity interpretation, despite the lack of unitary dynamics on 'Hphys- Another example of a 
self-adjoint operator on 'Hphys ls provided by the time-dependent Hamiltonian, which we 
shall study next. 



6. Self-adjointness of the Hamiltonian(s) 

Here it is shown that each element of the 1-parameter family of Hamiltonians for 
the classical Gowdy model can be promoted to a self-adjoint operator on TL and/or on 
Hphyg. At first sight, this result seems to contradict the failure of time evolution to be 
unitarily implemented. There is no contradiction, however. The usual link between unitary 
transformations and self-adjoint generators (Stone's theorem [26]) is, precisely, that every 
continuous 1-parameter unitary group has a self-adjoint generator, and vice versa. Because 
the classical Hamiltonian depends explicitly upon the time t, the dynamical evolution 
(2.23)-(2.28) does not form a 1-parameter group and Stone's theorem does not apply. 

The classical Hamiltonian for the polarized Gowdy model can be viewed as a 1- 
parameter family of functions H(t):T — > R, which can be expressed in any of the following 
equivalent forms: 

H(t) = lE(t) 

= Yt^ + ^ \ 2 Afl ^ ana - n + B n(t)^n a n + ^A* n (t)a* n a*_ n 
where 

A n (t) = ^n 2 [(H (\n\t)) 2 + (#i(|n|t)) 2 ] (6.2) 

Bn(t) = f n 2 [\H (\n\t)\ 2 + |#i(|n|t)| 2 ] . (6.3) 

The functions H(t) are generators of the transformation (2.23)-(2.28) in the following 
sense. The infinitesimal form of the transformation (2.23)-(2.28) at time t is given byf 

SI = = <hu SQ = = 8P, (6.4) 

d 

S<f = <P,t, S(t<P,t) = ^W,t)- ( 6 - 5 ) 



f Note that the variation does not commute with the time derivative since the form of the 
canonical transformation (like its generating function) depends upon time. 
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We have, for any tangent vector (5l,5jo,5<p), 

u,(6l,6'yo,6 ( p;0,0,<P,t) = SH(t). (6.6) 



To define H(t) as a family of operators on H,, we proceed as we did when defining the 
constraint operator C. For each fixed time t = r we compute the 1-parameter symplectic 
group on (T, Co) generated by H(t).* As we shall see, this one parameter symplectic group 
can be implemented as a continuous unitary group on Ti. The infinitesimal generator - the 
Hamiltonian - can then be defined via Stone's theorem. Finally, it is easy to check that 
this unitary group preserves 'H p hy S ^ so that the Hamiltonians thus defined are self-adjoint 
operators on 'Hphys- Here are the details. 

The 1-parameter group generated by H(r) can be viewed as a transformation on T 

(Z,70»9>P»Ojfc»Ojfe) — ► (Ks),lo(s),q(s),p{s),a k (s),al(s)) (6.7) 
defined by ("'" = £) 

Cd(8l, S-fo, <fy>, 1, 70, (p) = SH(t), (6.8) 

so that 

Q(s) = P(s) = p(s) = 0, ?(*) = — (6.9) 



a m (s) = -iB m (r)am(s) - iA* m (T)a*_ m (s) 
a* m (s) = iB m {T)a* m {s) + iA m (r)a- m (s). 



(6.10) 
(6.11) 



The point particle degrees of freedom (Q, P) are group invariants. The zero- frequency field 
modes (q,p) transform as do the coordinate and momentum of a free particle with mass r 
under time evolution. These transformations are certainly implementable as a continuous 
one-parameter unitary group. The transformation of the non-zero frequency field modes 
remains to be considered. The solution of (6.10)-(6.11) is given by 



a m (s) = 
fl m(«) = (%(«))*■ 



/i i \ -BrnKT) . n I \ 

cos(|m|s) — i — ; — ; — sm(|m|s) 



m 



Om(0) 



a ™ (T) sin(|m| S )a*_ m (0), 



m 



(6.12) 



Using the theory of unitary implementability [24,22] of symplectic transformations, the 
transformation (6.12) is, for each s, unitarily implementable if and only if 



E 



A m (r) . 



m 



sin(|m|s) 



< oo. 



(6.13) 



We emphasize that this group is not the set of time evolution canonical transformations. 
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It is straightforward to check that 

l^m(r)| 2 = ^2+0(^), (6.14) 

so that (6.13) is satisfied. 

Furthermore, the transformation (6.12) can be implemented as a continuous, unitary, 
1-parameter group if it is strongly continuous in the norm (3.14) [15]. Strong continuity 
means 

lim \\a(s) - a(so)\\ = lim > |a n (s) - a n (s )| =0, (6.15) 

which is easily verified as follows. The Bogoliubov coefficients in (6.12) are bounded so 
that there is an n and s-independent constant such that 

|on(s)| < (const.) (\a n (so)\ + \a- n (s )\). (6.16) 

Therefore 

\a n (s) - a n (s )\ 2 < (\a n (s)\ + |a n (s )|) 2 

2 2 (6-1'j 

< (const. )[\a n (s )\ + |a_ n (s )| + |a_ n (s ) | |a n (s ) |]- 

The right hand side of this inequality defines a square summable sequence of real numbers, 
thanks to the square-summability of a n (so). (Square summability of the first two terms is 
obvious, the last follows from the Schwarz inequality.) By the Weierstrass M-test [27] this 
guarantees that Yl \ a n( s ) — a n (so)| 2 converges uniformly for all s. Uniform convergence 
guarantees that \a n (s) — dn(so)| 2 converges to a continuous function of s, implying 
(6.15). 

From all these considerations, (6.12) is implementable as a continuous unitary group 
U(s):H^H, 



u\s)a m U(s) = 



n in, -Bmyr) ... . , 
cos(|m|sj — i — : — : — sm(|m|sj 
m 



i m --^P-sm(\m\s)al m , (6.18) 



from which the Hamiltonian H(r) is uniquely defined, up to an additive multiple of the 
identity, as the infinitesimal generator. 

The continuous unitary group U(s):H — > H acts on the physical Hilbert space.* We 
can then say that the Hamiltonian H(t) represents an observable (although it is not known 
precisely what is the domain of H(r) in Hphyg). To see that U(s) acts on H p h ys: we first 



* Formally, this follows from the fact that C and H(r) commute. But without a precise 
characterization of the domain of H(r) it is hard to conclude anything from this formal 
result. 
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note that the image of the state defined in (4.5) is a physical state, U(s)^q e "Hphys- 
Indeed, writing (6.18) as 

U\s)a m U{s) = a m (s)a m + /? m (s)al m . (6.19) 
a straightforward computation shows that 

U(s)* = N(s) exp J -i^-s ln{s)a{al n \ tf , (6.20) 

where 

a n (-s) 

and N(s) is fixed (up to a phase) by normalization. Evidently, the action of U(s) on 
the Fock vacuum is given by "pair creation" with each pair having zero total momentum, 
thus yielding a state in Hp^yg. Using (6.19) and (6.20) it is straightforward to compute 
the action of U(s) on a vector obtained as the image of any polynomial in the creation 
operators applied to ^q. Since these states span Ji (as x varies over a basis for L 2 (R 2 )), 
this defines U(s). It is then easy to see that states satisfying (3.20) are mapped by U(s) 
into states satisfying (3.20). For example, 

U(s)ala f _ n ^o = a* 2 (-s)aL n o4 

t t 2 1 (6 - 21) 

+ a^(-s)p*(-s)(a l n a n + a- n a l _ n ) + (3* (-s)a_ n a n U(s)^ , 

from which it is clear that 

,t„t 



U(s)aW- n ^0 e T-iphys- 



7. Remarks on the physical viability of the model 

In the quantum mechanical description of systems with a finite number of degrees of 
freedom, lack of unitary dynamics is normally associated with a failure of the probability 
interpretation of the model. The absence of unitary time evolution in the Gowdy model 
is an ultraviolet effect of the same sort as observed in [16,17]; it has no analog in the 
quantum mechanics of a system with a finite number of degrees of freedom. Indeed, we 
have seen that, despite the lack of unitary dynamics, the probability interpretation of the 
quantum Gowdy model appears to be intact in the following sense. The basic dynamical 
variables (Q, P) and (</?, (pj) (with the latter smeared with smooth functions of x) are 
self-adjoint operators on H for all t > 0. Observables (self-adjoint operators on H-phys) - 
besides functions of (Q, P) - can be built from the field variables (<£>, <pj)t by projection. 
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The spectral theorem then guarantees that the set of possible outcomes of a measurement 
of the observables have probabilities which add up to unity for all t > 0. Remarkably, 
even more complicated observables such as the Hamiltonians H(t) for the model can, 
for each t > 0, be defined as self-adjoint operators and given a consistent probability 
interpretation. Similar remarks can be made about the systems considered in [16,17]. 
There, the basic linear fields can, for all time, be defined as self-adjoint operators, so that 
the usual probability interpretation is available for them, despite the fact that the time 
evolution is not unitarily implemented. 

In [28] it is argued that non-unitary Schrodinger picture evolution in quantum gravity 
leads to difficulties with causality and locality. It is also pointed out there that these 
difficulties may be absent in a formulation of dynamics in the Heisenberg picture. Of course, 
in quantum mechanics the mathematically distinct Schrodinger and Heisenberg pictures 
are physically equivalent. But this is precisely because of the unitary implementation 
of dynamics in either picture. In the type of situation being discussed in this article, 
dynamical evolution is defined in the Heisenberg picture by the field equations, and the 
Schrodinger picture description of dynamics is unavailable. In this way the Gowdy model, 
as well as the models appearing in [16,17], appear to evade the unacceptable behavior 
discussed in [28]. 

For these reasons, it appears that the lack of unitary dynamics need not render the 
quantum Gowdy model physically unacceptable. Still, an important question left open 
by these considerations is the status of fundamental geometrical quantities in the Gowdy 
model, e.g., the metric and curvature . Some of the features of the metric and curvature 
operators have been studied in [3,4,10], but a more thorough investigation is warranted. To 
give a flavor of the issues involved, let us consider the metric components in (2.1) from the 
point of view of the quantum theory. Evidently, the exponentials of the quantum fields ip 
and 7 are required to define the quantum metric. Of course, it is too much to ask that the 
quantum metric components be defined by self-adjoint operators pointwise; some smeared 
version is required. Since Q and <f)(f) (defined in (5.3)) are self-adjoint, we can define a 
smeared self-adjoint <p via 

¥>(/) = e-sfyU), (7.1) 

which can be exponentiated to define a self-adjoint, smeared metric component. To define 
the quantum g yy and g zz metric components in an arbitrarily small neighborhood of any 
point one need only choose a sufficiently well-localized smearing function /. The definition 
of the metric function 7 (needed for gtt and g X x) is more problematic. The field 7 is 
formally defined in terms of (Q, P, 0, P$) via eqs. (2.5), (2.12), (2.13), and (2.20). It is 
not clear that smearing 7 with a smooth function of x and, e.g., normal-ordering of the 
creation and annihilation operators, will be adequate to render the expressions involving £ 
and n well-defined. That question aside, it is not at all clear that the variable 70, defined 
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in terms of (Q, P,(f>, P^) via (2.13) is well-defined. For example, the integral appearing in 
the last term of (2.13) is the generator for the following 1-parameter group of canonical 
transformations 

0( a ) = ^ P ^ p^ a ) = e -«p f (7.2) 

It is straightforward to verify that this transformation group is not unitarily implemented 
on H, precluding the existence of a self-adjoint generator. Evidently if 7 can be defined 
at all it will be through some sort of regularization procedure yet to be constructed (cf. 
[10]). It would seem that the need for a physically tenable definition of 7 in the quantum 
theory is the salient difficulty with the Gowdy model, not the failure of unitarity. 

Supposing one could not satisfactorily resolve the question of how to define 7 in the 
quantization of the Gowdy spacetimes under consideration, one might respond by searching 
for a different representation of the canonical commutation relations (CCR) than that used 
thus far, one which allows for a well-defined metric operator. Moreover, it is possible that 
other representations of the CCR will also allow for unitarily implemented time evolution.* 
However, the need for a search for an optimal choice of representation can be eliminated 
if one is willing to define the quantum Gowdy model using the algebraic quantum field 
theory (AQFT) formalism, such as been advocated by Wald [22], who was motivated by 
issues arising in quantum field theory in curved spacetime. In this formalism one, in 
effect, uses all representations of the CCR at once. Questions of unitary implementability 
become moot in the algebraic approach since they are tied to a particular representation, 
which is physically irrelevant from the point of view of AQFT. Indeed, it is sufficient to 
demonstrate that time evolution defines an automorphism of the C* algebraic structure 
used to define the theory, in which case the observables appearing in the C* algebra will 
be defined as self-adjoint operators for all time in any representation. An application of 
AQFT to the Gowdy model will have to take account of the non-linear nature of the space 
of solutions to the field equations - the constraint in particular, and the need to give 
probability distributions for the metric, curvature, etc. This certainly seems feasible and 
will be pursued elsewhere. 
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* This was suggested in reference [11]. Some of Berger's work [3] can be interpreted as 
considering alternate representations of the CCR. 
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